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Abstract: We introduce a new jet algorithm called XCone, for exclusive Cone, which is 
based on minimizing the event shape Y-jettiness. Because Y-jettiness partitions every event 
into Y jet regions and a beam region, XCone is an exclusive jet algorithm that always returns 
a fixed number of jets. We use a new “conical geometric” measure for which well-separated 
jets are bounded by circles of radius R in the rapidity-azimuth plane, while overlapping jet 
regions automatically form nearest-neighbor “clover jets”. This avoids the split/merge criteria 
needed in inclusive cone algorithms. A key feature of XCone is that it smoothly transitions 
between the resolved regime where the Y signal jets of interest are well separated and the 
boosted regime where they overlap. The returned value of Y-jettiness also provides a quality 
criterion of how Y-jet-like the event looks. We also discuss the Y-jettiness factorization 
theorems that occur for various jet measures, which can be used to compute the associated 
exclusive Y-jet cross sections. In a companion paper [1], the physics potential of XCone is 
demonstrated using the examples of dijet resonances, Higgs decays to bottom quarks, and 
all-hadronic top pairs. 
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1 Introduction 

Collisions at the Large Hadron Collider (LHC) are dominated by jets, collimated sprays of 
hadrons arising from the fragmentation of energetic quarks and gluons. Jets are crucial to 
connect the observed hadronic final state to the short-distance hard interaction. Fundamen¬ 
tally, the definition of a hadronic jet is ambiguous, since there is no unique way to map 
color-singlet hadrons to color-carrying partons. Moreover, different physics applications can 
benefit from different jet definitions. For these reasons, a wide variety of jet algorithms have 
been proposed to identify and study jets [2, 3], though currently, most LHC measurements 
involve jets clustered with the anti-Ur algorithm [4], 
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In this paper, we present a new jet algorithm that we call “XCone”. It is based on 
minimizing the event shape IV-jettiness [5] and uses developments from the jet shape N- 
subjettiness [ 6 , 7]. The key feature is that IV-jettiness defines an exclusive cone jet algorithm. 
Like the exclusive for algorithm [ 8 ], our XCone algorithm returns a fixed number of jets, 
relevant for physics applications where the number of jets is known in advance. Like anti-for 
jets [4], XCone jets are nearly conical for well-separated jets, such that they have fixed active 
jet areas [9, 10]. Typically, when using other jet algorithms, the boosted regime of overlapping 
jets requires separate analysis strategies using fat jets with substructure [11-14], In contrast, 
with XCone the jets remain resolved even when jets are overlapping in the boosted regime. In 
this way, XCone smoothly interpolates between the resolved regime of widely-separated jets 
and the boosted regime of collimated subjets. This feature will be explored in more depth 
in a companion paper [ 1 ], which demonstrates the application of XCone for the examples of 
dijet resonances, Higgs decays to bottom quarks, and all-hadronic top pairs. 

The possibility of using IV-jettiness as a jet algorithm was already pointed out in ref. [5] 
and further explored in ref. [7]. Here, we more fully develop the idea of IV-jettiness jets 
and present a concrete implementation of the XCone algorithm. As a global event shape, 
IV-jettiness measures the degree to which the hadrons in the final state are aligned along 
N jet axes or the beam direction. It was originally introduced to veto additional jets in an 
event, providing a way to define and resum exclusive IV-jet cross sections [5, 15, 16]. 1 IV- 
jettiness was later adapted to the jet shape IV-subjettiness [ 6 ], which is an efficient measure to 
identify IV-prong boosted hadronic objects such as top quarks, W/Z bosons, and Higgs bosons 
within a larger jet (see also [17]). By minimizing IV-(sub)jettiness, one can directly identify 
N (sub)jet directions, and a fast algorithm to perform this minimization was presented in 
ref. [7]. IV-jettiness jets have been used to resum the invariant mass of nearby jets [18], to 
make predictions for jet mass spectra [19, 20], for studying DIS and nuclear dynamics [21-26], 
and to define recoil-free jet observables [27]. As an IV-jet resolution variable, IV-jettiness has 
been utilized to combine perturbative calculations with parton showers in Geneva [28], and 
very recently to define a powerful subtraction scheme for fixed-order calculations at next-to- 
next-to-leading order [29, 30]. 

As we will see, there is considerable flexibility in precisely how one defines IV-jettiness, 
and several different IV-jettiness measures yielding different jet regions have been considered 
before [5-7, 16, 19]. Here, as the XCone default, we propose a “conical geometric” measure 
that incorporates the insights from the different previous use cases. This measure is based on 
the dot product between particles and lightlike axes as in ref. [5] but incorporates an angular 
exponent f3 as in ref. [7], as well as a beam exponent 7 for additional flexibility (see table 1 

1 The reader should be aware that there are two different definitions of “exclusive” which are both standard 
in their respective contexts. An exclusive IV-jet algorithm is one that returns exactly N jets, regardless of 
what happens in the rest of the event. An exclusive IV-jet cross section is the rate to produce exactly N jets, 
with a restriction on what happens in the rest of the event. XCone is an exclusive IV-jet algorithm, but it 
can be used either to measure inclusive IV-jet cross sections (if there are no restrictions made on unclustered 
particles) or an exclusive IV-jet cross section (if there is a restriction, say, that 7Iv < 7f u t). 
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below). Crucially for the purposes of jet finding at the LHC, this measure yields (nearly) 
conical jets over a wide rapidity range, and the user can choose the desired jet radius R. 

For most physics applications, we propose a default setting of (5 = 2 and 7 = 1 , which 
acts similarly to existing cone algorithms (see e.g. [31-34]) in that the resulting jet regions 
are (approximately) stable cones where the jet momenta and the jet axes align. The key 
difference to algorithms like SISCone [34] is that XCone does not require a split/merge step. 
In particular, typical inclusive cone algorithms have an overlap parameter which determines 
whether two abutting stable cones should be joined or remain separate. By contrast, XCone 
only requires setting the jet radius R and the number of desired jets N, and the split/nrerge 
decision is determined dynamically through IV-jettiness minimization. In a companion paper 
[1], we show examples of quasi-boosted kinematics that capitalize on this exclusive approach 
to cone jet finding. 

There are interesting connections between IV-jettiness minimization and previous work 
to define jets via cluster optimization [33, 35-43]. Stable cone finding is closely related to 
1-jettiness minimization with (5 = 2 [33], and similar algorithms are relevant for a recently 
proposed “jet function” 2 3 optimization strategy [47-49]. One can even prove an exact equiv¬ 
alence between these algorithms when finding a single cone jet of fixed opening angle [50]. 
Finding the thrust axis [51] is related to 2-jettiness minimization with (5 = 2. 3 There is also 
an observable called triplicity [53] which is related to 3-jettiness. For a general N, fc-means 
clustering [54] (with k = N ) is a type of IV-jettiness minimization, with (5 = 2 correspond¬ 
ing to traditional fc-means and (5 = 1 corresponding to Rl-fc-means [55]. In all these cases, 
IV-jettiness minimization is an infrared and collinear (IRC) safe procedure. 

Because cluster optimization is a difficult computational problem, our practical XCone 
implementation will use recursive clustering algorithms [8, 56-59] to approximate IV-jettiness 
minima. Roughly speaking, we run a generalized kx clustering algorithm to determine IRC- 
safe seed jet axes as a starting point for an iterative one-pass minimization algorithm, in which 
IV-jettiness is used to find the final jet axes and define the jet regions. Separating jet axes 
finding from jet region finding appeared previously in the context of recoil-free jets [27, 60], 
where a fixed radius cone was centered on winner-take-all axes [27, 61, 62] or broadening 
axes [7, 27]. XCone allows us to extend this strategy to IV-jet events, with (5 = 1 yielding 
recoil-free jets and (5 = 2 yielding traditional cones where the jet axes and jet momenta are 
(nearly) aligned. 

A key feature of the measures we consider, including the default XCone measure, is that 
IV-jettiness can be decomposed into a direct sum of contributions from the jet and beam 
regions. When utilizing measures with this property, there exist active-parton factorization 
theorems for IV-jettiness cross sections valid to all orders in a s . Furthermore, the default 

2 The name jet function in this context should not be confused with the more standard usage in the context 
of factorization of cross sections into hard, soft, and jet functions, e.g. [44-46]. Here our primary use of the 
name jet function will be in this factorization context, see sec. 5. 

3 Naively, one might think that spherocity [52] should be related to 2-jettiness with /3 = 1. However, 
minimizing this quantity does not give rise to the spherocity axis, but rather to kinked broadening axes [27]. 


- 3 - 



XCone measure is linear in the particle momenta which greatly simplifies the calculation of the 
perturbative jet and soft functions needed to determine the IV-jettiness cross section. Thus, 
the ingredients needed for higher-order logarithmic resummation or fixed-order calculations 
are simpler for jets defined with the XCone algorithm, in contrast for example to those defined 
with clustering algorithms like anti -Jct- We will discuss these factorization theorems in some 
detail for various choices of IV-jettiness measures, including the XCone default. 

The remainder of this paper is organized as follows. In section 2, we review how to define 
an exclusive jet algorithm via minimizing IV-jettiness. We then discuss a variety of IV-jettiness 
measures in section 3, including the conical geometric measure that is the basis for XCone. 
In section 4, we discuss some details of our XCone implementation, in particular the choice of 
seed axes for finding a (local) IV-jettiness minimum. In section 5, we discuss the factorization 
theorems for IV-jettiness with various measures. This section is more theoretically technical 
than the others and may be skipped by readers not interested in this factorization. We con¬ 
clude in section 6. The XCone algorithm is available through the Nsubjettiness Fast Jet 
Contrib [63, 64] as of version 2.2.0. 

2 IV-jettiness as a Jet Algorithm 

Given a set of normalized lightlike axes ha = {1 ,hU} with n\ = 1, IV-jettiness is defined as 4 

f N = ^mm{pj et {p i ,ni),...,pj et {p i ,n N ),p heain (pi)} . (2.1) 

i 

The sum runs over the four-momenta pi of all particles that are considered as part of the 
hadronic final state and should take part in the jet clustering. The Pj e t(Pi,nA) is a distance 
measure to the A-th axis ua, and Pbeam(Pi) is a distance measure to the beam. Depending on 
the context, the beam measure can be separated into two beam regions with lightlike beam 
axes n a ,b and (partonic) center-of-mass rapidity Y such that 

Pbeam (jPi ) => min{pbeam(Pi, n a , Y ), p b eam(Pi, n b , T)} . (2.2) 

This form will be relevant for the discussion in section 5. 

For a given form of pj e t and Pbeami the minimum inside Tn in eq. (2.1) partitions the 
particles i into N jet regions and an unclustered beam region. To use IV-jettiness as a jet 
algorithm, one minimizes Tn over all possible lightlike axes directions: 

Tn = min Tn- (2.3) 

ni,7i2,...,rijv 

The locations of the axes at the minimum define the centers of the jet regions. In previous 
applications, one uses a separate method to choose the IV-jettiness axes ha, e.g. from the N 
hardest jets found by some other jet algorithm. One then uses Tn only for the jet partitioning 
(in which case there is no need to distinguish Tn = Tn)- This use of Tn already provides a 

4 Here we use a dimension-one definition as in refs. [16, 19] instead of the dimensionless tn used in ref. [5]. 
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well-defined and IRC-safe way to define N exclusive jets. The additional overall minimization 
in eq. (2.3) over the axes ua promotes 7/v to a standalone exclusive jet algorithm. This axis 
minimization is nontrivial and we discuss our strategy to perform it in section 4. 5 Note 
that “minimization” can refer either to finding the global 7/v minimum or using an IRC-safe 
procedure to find a local 7/v minimum, either of which is suitable for the discussion below. 

Any choice of measure together with the specific algorithm to minimize 7/v defines an 
exclusive jet algorithm. In particular, 7/v in eq. (2.1) always identifies N jet regions (and 
one beam region), regardless of how close the axes ua might be to each other. When the 
axes are well separated, the boundary of the jet regions is determined through competition 
between pj et and Pbeam- When the axes are close together, the jet regions are determined by 
the competition between different pj et . 

To go from an exclusive jet algorithm to an exclusive cone jet algorithm (i.e. XCone), one 
wants the jet boundaries to approximate circles in the rapidity-azimuth plane, which can be 
achieved by an appropriate choice of jet and beam measures. In section 3, we study a variety 
of jet and beam measures which are summarized in table 1. This includes three new measures: 
the general conical measure in eq. (3.3) which yields exact cones for widely-separated jets; 
the modified geometric measure in eq. (3.10) whose jet measure is linear in particle momenta 
like the original geometric measure but exhibits smooth behavior at zero rapidity; and the 
recommended XCone default in eq. (3.18) which yields approximate cones and also features 
this linearity. By construction, the XCone default measure yield jets with approximately 
fixed active jet areas over a wide range of jet rapidities. 

In addition to partitioning the event into jet and beam regions, the returned value of 
T/v is a quality criterion that measures how well an event is characterized by N jets. The 
contribution to the 7/v value from a given jet provides a measure of how collimated the jet is. 
For narrow jets (i.e. small effective jet radius), 7/v is typically dominated by the contribution 
from the beam region. Thus, for LHC applications, one typically wants pbeam(Pi) to be 
proportional to pTi (the transverse momentum of particle i ) such that minimizing 7/v results 
in the least unclustered px■ Larger values of 7/v, and its beam contribution in particular, then 
indicate additional activity or hard jets in the event. An improved measure of jet quality can 
be obtained by examining the individual jet and beam contributions to 7/v, as in [16]: 

N 

r A r = r^ eam + r^ ets = r^ eam + ^r^. (2.4) 

A= 1 

Here, 7y tb provides a global measure for assessing how collimated the jets are without con¬ 
tamination from the beam region, and one can obtain individual quality measures for each of 
the N jets by examining their individual numerical contributions 7"[y to the total X-jettiness. 
In section 5, we discuss some of the theoretical aspects involved in calculating 7/v as well as 
the cross section that is fully differential in Ty eam and the N observables 7jy. 

5 One might also be able to dynamically determine the total rapidity Y or the beam axes n a: b through 
minimization, though that feature is currently not present in the XCone code. 
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Name 

Pjet ijPii T^Al) 

Pbeam (Pi) 

A « ttR 2 1 

Conical [ 7 ] 

(RiA\P 

PTi 

PTi 

/ 

General Conical 

pnf{pi ) y^-y 

PTi f(Pi ) 

/ 

Geometric [ 19 ] 

UA ' Pi 

Po 

TOTie -iM 


Modified Geometric 

UA 'Pi 

mn 


Po 

2 cosh iji 


Geometric-i? [ 19 ] 

a A ■ Pi 

p{R,Va) 

rriTie~\ Vi \ 

/ 

Modified Geometric -R 

UA - Pi 

m T i 

/ 

pc{RiVa) 

2 cosh 7 /i 

Conical Geometric 

PTi / 271 A -Pi 1 \P /2 

PTi 

/ 

( 2 cosh UiT~ l \UTAPTi R 2 ' 

(2 coshT/i ) 7-1 

XCone Default (/3 = 2,7 = 1 ) 

2 cosh U a 

R2 UA-Pi 

PTi 

/ 

Recoil-Free Default (/? = 1,7 = 1 ) 

1 2 cosh yA 

Y R 2 P Ti nA • ' (H 

PTi 

/ 

II 

1 ° 

II 

to 

cosh 7 /A 

PTi 

/ 

cosh 7/i I? 2 77 ’ ? 

2 cosh ?/,: 


Table 1: IV-jettiness measures studied in this paper. The conical geometric measure with 
/3 = 2 and 7 = 1 is the suggested XCone default, giving stable cone jets (like the conical 
measure) through dot-product distances linear in pi (like the geometric measures). The recoil- 
free variant with /? = 1 centers the jet around its hardest cluster, making the jet regions less 
sensitive to soft contamination. In the conical geometric measure, rip a = 1/cosh y\. In 
the (modified) geometric -R measures, p(c){R,yA) is a rapidity-dependent scale factor that 
yields jet areas of exactly irR 2 (though not conical jet boundaries). The checkmarks indicate 
measures that yield jets with active areas of ~ ttR 2 for well-separated jets. These active areas 
are ttR, 2 to within < 1% over a wide rapidity range (see figure 4 below). 


Before discussing the specific measures, we want to make a general comment about un¬ 
derlying event and pileup, two effects that are known to impact jet reconstruction. While the 
value of Tn depends strongly on these effects, the jet regions found by minimizing Tn are 
no more sensitive to underlying event and pileup than traditional jet algorithms. The reason 
for this mismatch is that the beam contribution to the Tn value can get large contributions 
from these effects, but the change in Tn as the axes ua are varied only depends on hadrons 
in the vicinity of the jet regions. This is particularly true for recoil-free measures, where 
the minimized axis direction is almost entirely insensitive to soft contamination [60]. For 
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pileup specifically, the minimization in eq. (2.3) remains sensible even with negative energy 
particles, so one has the option of introducing negative energy ghosts as a way to implement 
area subtraction [9, 10, 65]. For isolated jets, one can derive a closed-form integral expression 
for the active jet area, which depends only mildly on the jet rapidity. 


3 Choice of Measure 


As already mentioned, every choice of jet and beam measure defines some kind of TV-jettiness 
jet algorithm. We now review previous measures in the literature en route to explaining 
the logic behind the new XCone default measure. Example jet regions found from some of 
these measures are shown in figures 1 and 2 for a boosted top event from the BOOST 2010 
event sample [11]. In figure 3, we show a comparison between the XCone default and the 
anti -kx algorithm [4] . While XCone and anti-fop are very similar for widely separated jets as 
in figure 3b, they behave quite differently when the jets are close together as in figure 3a. A 
more extensive discussion and anti-fcp comparison can be found in the companion paper [1], 


3.1 The Conical Measure 

The first conical IV-jettiness measure was proposed in ref. [7]: 6 


Conical Measure 


Pjet ('Pi ; Tf ) VTi 



Pbeam (Pi ) — VTi j 


(3.1) 


where 

RiA = V(Vi - Va) 2 + (<fc - 4>a) 2 (3.2) 

is the distance between p, and ha in the rapidity-azimuth plane, and /3 is an angular weighting 
exponent. The parameter R acts like the jet radius in a cone algorithm, since particle i can 
only be clustered into jet A if p^t{Vii n A) < Pbeam(Pi)) which is equivalent to Ria < R■ Thus, 
the measure in eq. (3.1) yields jets that are exact circles with radius R in the rapidity-azimuth 
plane, as shown in figure 2a, unless two jet axes are closer than R. When two or more axes are 
closer than R to each other, the jet regions are determined by Voronoi partitioning (i.e. nearest 
neighbor). This yields “clover jet” configurations as shown in figure la. 

For small R, Tn is dominated by the beam measure, which is just the unclustered pr in 
an event. Thus, this measure typically finds the N jets with the largest vt in an event. By 
adjusting the exponent /3, the jet axis can be varied to point along the jet direction (/3 = 2, 
“mean”) or along the hardest cluster inside a jet (/3 = 1, “median”), see also refs. [7, 27, 60]. 

Naively, the conical measure might seem to be the only measure yielding conical jets, 
since any change to the measure would affect the competition between pj et and Pbeam and 

6 Strictly speaking, the measure in ref. [7] has an extra rapidity cut parameter. 
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Geometric Measures 


Conical Measure (J3 = 2) 



y y 

(a) (b) 


XCone Default (Conical Geometric) 



(c) 


Figure 1: Jet regions found with various iV-jettiness measures. This is a tt event from the 
BOOST 2010 event sample [11], and every measure has N = 6 and R = 0.5. (a) Conical 
measure with (5 = 2. (b) Original and modified geometric measures, (c) Conical geometric 
measure with (5 = 2 (XCone default) and (5=1 (recoil-free default). The conical and conical 
geometric measures yield (approximately) circular jets. For all measures, the overlap region 
between jets is automatically partitioned by nearest neighbor, as given by the jet measure. 





Conical Measure (J3 = 2) Geometric Measures 



XCone Default (Conical Geometric) 



(c) 


Figure 2: Same tt event as in figure 1, but for N = 2 and R = 1.0. (a) The conical measure 
yields exactly circular jet regions for widely-separated jets, (b) The geometric measure ex¬ 
hibits cusps at y = 0 which are smoothed out with the modified geometric measure, (c) The 
XCone default (/? = 2) yields jets centered along the total jet momentum while the recoil-free 
default (/? = 1) yields jets centered along the hardest cluster within the jet. 
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Figure 3: Comparison between the XCone default (/3 = 2) and anti -fop, using the same tt 
events as figures 1 and 2. (a) Unlike anti-fop which merges jet regions closer in angle than 
« R, XCone allows such jet regions to remain split, (b) For widely-separated jets, XCone 
yields nearly identical jet regions to anti-fop. 


change the style of the event partitioning. One can maintain conical jets, however, if one 
deforms eq. (3.1) via 


General Conical Measure 


Piet{pi,n A ) =PTif(Pi) 



Pbeam (.Pi) PTi f (jPi) i 


(3.3) 


where f(pi) is any dimensionless function of the particle four-momentum. This measure still 
returns exactly conical jets with overlapping jets still having Voronoi partitioning, because 
the factor of f(pi) drops out when comparing p ]et to Pbeam or when comparing two different 
Pjet- While the partitioning for given axes does not depend on f(pi), the f(pi) factor does 
play a role in determining the overall Tn minimum in eq. (2.3). So the final jets will have 
different axes depending on the choice of f(pi). We will exploit this possibility when defining 
the conical geometric measure in section 3.3. 


3.2 The Geometric Measure 


A variety of IV-jettiness measures were proposed and studied in refs. [16, 19]. For the purposes 
of defining a cone jet algorithm, the most promising choice is the geometric measure: 


Pjet (Pi ; a) 
Pbeam (Pi) 


nA •Pi 
Po 

min{n a • Pi, n b ■ p t } , 


Geometric Measure 


(3.4) 









where n a ,b = {1, 0,0, ±1} and the z-direction is the beam direction, such that 

min {n a • p h n b ■ p^ = p° - \p?\ = m T ie~ M . (3.5) 

Here, mTi = \JPj^ + m 1 , Vi is the rapidity, and this is the form given in table 1. 

The presence of the n-pi dot product in the jet and beam measures is very natural from a 
theoretical perspective, since it makes the measure linear in both pi and n. The linearity in the 
jet axes n A implies that the total jet three-momentum is exactly aligned with the axis direction 
n A (see section 4.2). The linearity in pi implies simple factorization properties for 7/v and also 
tends to make perturbative calculations much simpler (see e.g. refs. [16, 18, 19, 25, 30, 66]). 
For this reason all IV-jettiness calculations so far which involve initial state hadrons have been 
based on measures linear in pi, like the geometric measure. 

Despite the presence of the dot product ha ■ Pi, the geometric measure actually behaves 
quite similarly to the conical measure.' To see this, note that the momenta pi and lightlike 
axes ha can be expressed as 


Pi = {m Ti cosh yi, p Ti , m T i sinh } , 
ua= { 1, n T A, tanhy A } , 

and their dot product is given by 


PTi = \pn\ , 

n TA = | nTA | = 7- , 

cosh y A 


nA ■ Pi rnn 


cosh (yi - i/a) - cos (<fo - (f> A ) 

uta PTi PTi 

In the limit of small angles and for massless particles we thus have 


(3.6) 

(3.7) 


(3.8) 


n A ■ Pi PTi R 2 iA 

, n A ) = -~ x-r-— • 

p 0 2 cosh y A Po 

Hence, the /3j e t for the geometric measure acts similarly to the general conical measure in 
eq. (3.3) with /3 = 2 and f(pi) = 1/(2cosh?/,), at least to the extent that coshy* ~ coshy A . 
This also shows that the parameter po in the geometric measures controls the size of the jet 
regions with roughly po ~ R 2 . 

Since the geometric measure does not take the precise form of eq. (3.3), it yields football¬ 
like jets in the central region with cusps at y = 0, which get accentuated for larger jets as 
shown by the green thick lines in figures lb and 2b. For overlapping jets, it produces similar 
clover jets due to the competition between the pj e t f° r different jets. 

Although not as extreme as the jet shapes obtained with an invariant mass measure 
(see ref. [16]), these cusps in the jet boundaries are somewhat unnatural for experimental 
applications. Since the shape of the jet regions is determined by the competition between /jj e t 

7 In the context of recursive clustering algorithms, this dot-product form was also mentioned as an option 
in ref. [8]. 



- 11 - 










and pbeam; we can modify the geometric measure to yield more conical jets by introducing an 
explicit compensating factor of f(pi) = 1/(2 cosh yj) in the beam measure: 


Modified Geometric Measure 


Pjet (Pi j 71 a) 


Pbeam (Pi ) 


nA •Pi 

Po 

m T j 

2 cosh i/i 


(3.10) 


With the approximations in eq. (3.9) and coshyj ~ coshyj this modified measure is now 
approximately the same as eq. (3.3) with ,3 = 2 and f(j>i) = 1/(2cosh yj). Hence, it yields 
reasonably conical jets also in the central region, as shown by the purple lines in figures lb 
and 2b. This corresponds to only a slight modification of the geometric beam measure, since 
close to the beam axes, i.e. for large yj, we have 

mT * - >m T ie~ M . (3. 

2 cosh yi 


This implies that the modified geometric measure has very similar factorization properties 
as the geometric measure, which we will return to in section 5. The use of 1/(2 cosh yj) to 
replace e - ^ is the same as the well-known distinction between using C'-parameter [67, 68] and 
thrust [51] event shapes to describe the narrow dijet limit in e + e - collisions, see e.g. [69-72]. 

While we can roughly associate p$ ~ R?, the jet area itself still differs from irR 2 , especially 
for larger R and away from central jet rapidities. To enforce jets of a constant jet area, 
regardless of the jet rapidity and jet boundary, ref. [19] also introduced a geometric -R measure 
where the jet measure is rescaled by a rapidity-dependent factor to maintain 7tR 2 jet areas 
for widely-separated jets: 


Geometric- R Measure 


Pi ‘ tl - P, ' nA] - PiR.M) 

Pbeam(.Pi) TTlTi6 ^ 


nA-Pi, 


(3.12) 


Here, p(R,ua) is given in terms of the the integral Io(a,/3 ) from [16] which determines the 
geometric jet area (for nonoverlapping jets) via the transcendental equation [20] 

f) += R2 • “± = 1± ‘™ h ^- < 313 > 

Numerical results for p were given in ref. [19]. The same modifications as above lead to the 
modified geometric -R measure 


Modified Geometric -R Measure 


Pjet (Pi j 7lA ) 
Pbeam ( Pi ) 


pc(R,va) nA ' Pi 

m T i 

2 cosh yj 


(3.14) 


where pc(R,Va) is different than in eq. (3.12) due to the difference in the beam measures. 

In all of the above cases, the rapidity suppression in the beam measures at large rapidi¬ 
ties makes Tn much less sensitive to the forward region. This means the Tn minimization 
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effectively corresponds to minimizing a rapidity-weighted sum of unclustered pt (he. the un¬ 
clustered beam thrust [15] or “beam C-parameter” contribution). As a result, the algorithm 
will dominantly identify central jets over forward jets, which could have interesting appli¬ 
cations, e.g. when one wants to avoid picking up forward jets from initial-state radiation. 
Corresponding forward-insensitive rapidity-weighted jet vetoes have been discussed recently 
in ref. [73]. 


3.3 The Conical Geometric Measure 


Combining the lessons of the conical and geometric measures, we now introduce the con¬ 
ical geometric measure which aims to combine their advantages. For a specific choice of 
parameters, this will be the XCone default measure. 

Like the conical measure, we want a measure that returns (nearly) conical jets, and we 
also want a parameter (3 in the jet measure to adjust the behavior of the jet axes. Like the 
geometric measure, we want a measure that depends on the dot products between lightlike 
axes and particles, since that is the simplest distance to use in theoretical calculations, and 
can be made linear in the particle momentum (here by choosing /3 = 2). These requirements 
lead us to 


Conical Geometric Measure 


Pjet (Pi j T A ) 
Pbeam (Pi ) 


PTi f %n A ■ Pi 1 
(2 cosh |/j) 7_1 V n T A PTi R 2 ) 
PTi 

(2 cosh t/j) 7-1 ’ 


where again uta = 1 / coshy^. In the jet measure, we recognize the last factor in parentheses 
as the approximate form for R{ A in eq. (3.8), which now yields jets that are very nearly 
conical. The (3 factor acts just like the (3 factor in the conical measure. For additional 
flexibility, we have chosen a common f(pi) = (2 cosh y*) 1-7 in the beam and jet measures. 
This multiplicative factor affects the axes found by minimization, but not the beam and jet 
regions. It is parametrized by 7, such that for 7 = 1 this reproduces the beam measure of the 
conical measure while for 7 = 2 this is closely analogous to the beam measure of the modified 
geometric measures. 

There is additional freedom in defining the conical geometric measure that we will not 
exploit in this paper. For example, we could multiply the jet or beam measures by any 
function of 

m T j 
PTi 

which would give slightly different behavior for massive hadrons. In the jet measure, we could 
multiply by any function of 

cosh IJi 
cosh y A ’ 

since this quantity is nearly one for narrow jets. For example, the modified geometric measure 
is reproduced exactly by taking f3 = 7 = 2 and in addition multiplying the beam and jet mea¬ 
sures by rriTi/pTi and coshy*/coshy^, respectively. These choices are somewhat analogous 


(3.17) 


(3.16) 
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Jet Area Comparison 



Figure 4: Comparison of the analytic jet areas for a single jet (N = 1). Unlike the modified 
geometric measure, the conical geometric measure (here shown for the XCone default of 
/3 = 2 ) has uniform jet areas as a function of rapidity. For R < 1.0, this area is within 1% of 
ttR 2 from the conical measure. 


to the choice of recombination schemes in recursive jet algorithms, since they are irrelevant 
for infinitely narrow cones and massless inputs. That said, for (3 = 2, 7 = 2 the factor of 
cosh da/ coshy, that appears in the conical geometric measure relative to the (modified) geo¬ 
metric measure ensures that the jet area is very close to irR 2 even for relatively forward jets, 
as shown in figure 4. 


3.4 The XCone Default Measure 


For LHC applications, our recommended XCone default is the conical geometric measure with 
(3 = 2 and 7 = 1 : 


XCone Default Measure ((3 = 2) 


Pjet (Pi > a) 


2 cosh ua 

^2 n A' Pi , 


PbeamiPi ) — PTi ■ 


(3.18) 


By choosing 7 = 1 , the beam measure is the same as the conical measure, so minimizing Tn 
minimizes the unclustered pr- By choosing (3 = 2. the jet axis (approximately) aligns with 
the total three-momentum of the jet, as is typical for traditional stable cone algorithms. Note 
that the jet measure is linear in pi , as desired for theoretical calculations. In figure 4 we show 
that the active area of XCone jets is very nearly ttR 2 for well-separated jets, see also ref. [1]. 

Alternatively, in cases where recoil-sensitivity [74-77] is an issue (such as in high pileup 
environments [60]) we can use (3=1 and 7 = 1 : 


Recoil-Free Default Measure ((3 = 1) 


Pjet (Pi 1 Ua) 


2 cosh yA 
R 2 


■ PTi UA- Pi : 


Pbeam(Pi) — PTi • 


(3.19) 
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Here, the jet center aligns approximately along the broadening axis of the jet [7, 27], which 
is the axis that minimizes the summed transverse momentum relative to it. This is similar 
to finding the “median” jet energy and the jet axis tends to point along the most energetic 
cluster within a given jet. Again, the jet area is approximately irR 2 . 

These XCone default measures are the basis for our LHC case studies in the companion 
paper [1], where we find that both f3 = 2 and (3 = 1 give comparable results for jet recon¬ 
struction (in the absence of jet contamination). The jet regions for XCone default are shown 
in figures lc and 2c. With a single energetic cluster inside a jet, the difference between (3 = 2 
and (3 = 1 is very small (again in the absence of jet contamination), analogous to the way that 
the mean and median of a peaked distribution are very similar. This is shown in figure lc. 
When a jet has substructure, the “mean” (/3 = 2) and “median” (/3 = 1) axes are offset, 
as shown in figure 2c for the same event with X = 2. One can also see that for larger jet 
radius, the jet regions are slightly elongated along the azimuthal direction compared to the 
rapidity direction. This arises because of the trigonometric functions in eq. (3.8). In ref. [1] 
it is mentioned that this deformation from exact circles yields slightly improved performance 
when reconstructing invariant-mass peaks. 

4 Details of the XCone Algorithm 

For a given X-jettiness measure entering in eq. (2.1), we need to implement the minimization 
procedure in eq. (2.3) to determine the jet axes ha- In general, the only guaranteed method 
to find the global minimum of Tn is to test by brute force all possible partitions of the 
final-state particles into X jet regions and one beam region. Since this is computationally 
prohibitive, our aim is to find good approximations of the global minimum by relying on 
methods that strictly speaking only find local minima of Tn. Even if the algorithm does not 
find a guaranteed global Tn minimum, as long as all steps are fully specified and IRC safe, it 
still represents a well-defined exclusive cone algorithm which retains the key features of the 
X-jettiness partitioning according to the specified jet and beam measures. 

Throughout this section, we restrict ourselves to the case 7=1, which is currently 
implemented in the XCone code and is also used by the default measures. 

4.1 One-Pass Minimization 

For the conical measure in eq. (3.1), ref. [7] introduced a modification of Lloyd’s method [54] 
that finds a local minimum of Tn for 1 < f3 < 3. We can adopt a similar strategy for more 
general measures. 

Our minimization algorithm proceeds as follows, with more details given below: 

1) Find seed axes: Determine a set of suitable IRC-safe initial axes ha- 

2) Assignment: For fixed axes ua , assign particles to jet and beam regions via Tn parti¬ 
tioning. 
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3) Update axes: For fixed partitioning, update axes ha via Tn minimization. 

4) If axes have converged then stop, otherwise go back to step 2). 

To be IRC safe, this procedure must be fully deterministic. We therefore always perform a 
one-pass minimization, i.e., the above algorithm is repeated precisely once per event without 
any stochastic elements (such as random variations in the seed axes). The procedure to 
determine the seed axes in step 1) is deterministic and IRC safe, as described in section 4.3. 
The seed axes are then iteratively improved to a local minimum of Tn in steps 2) and 3). 

In the assignment step 2), the final-state particles are assigned to one of the N jet regions 
or to the beam region via the Tn partitioning in eq. (2.1) for the current set of fixed trial 
axes n a • This step can be easily implemented for any choice of measure as it only depends 
on the competition between the jet measures p^tiTAiVi) for fixed ua and the beam measure 
Pbeam (pi)■, so we do not need to discuss it further. 

In the update step 3), the axes ua are improved to minimize the contribution to the Tn 
value within each jet region, keeping the jet constituents determined by the partitioning in 
the previous assignment step fixed. Different update steps are needed for different measures, 
since there is no general procedure to find the axes ua that minimize JT pj e t(riA, Pi) - 8 Once 
an appropriate update step is found, the assignment and update steps can be iterated until 
the axes converge to within some specified accuracy. In section 4.2, we describe a general 
update step that works well for the measures studied in this paper. 

As discussed in ref. [7], these one-pass minimization procedures are quite effective for 
IV-subjettiness, often converging to the global minimum. There are additional complications, 
however, for IV-jettiness. The reason is that IV-jettiness has a beam region, and particles in 
the bulk of the beam region are insensitive to small changes to the location of the jet axes 
ua- Even minimization routines that try to go “uphill” to escape local minima may never 
find the optimal jet axes. Given that Tn corresponds roughly to the unclustered pt in an 
event (for 7 = 1), failing to find a decent Tn minimum means that one will identify too many 
soft jets. Therefore, for XCone to be a practical jet algorithm, one has to find a good set of 
seed axes for one-pass minimization. In section 4.3, we show how to find such seed axes by 
utilizing recursive clustering algorithms. 

Another possibility to further improve the Tn minimization is by running the above (or 
any other) exclusive jet algorithm to find N + n jets. Starting from these, one can then 
perform the remaining partitioning into N jets by explicitly testing all possible combinatorial 
options to find the best minimum. This option is available in the XCone code, though not 
recommended by default for reasons of speed. One advantage of this strategy is that it reduces 
to the exact Tn minimization for up to N + n final-state particles. This makes it convenient 
for fixed-order calculations up to N”LO, where the one-pass minimization with seed axes 
could induce rather complicated boundaries in the phase-space integrations. 

8 Even if one does find such a procedure, one has to check on a case-by-case basis whether the assign¬ 
ment/update iteration actually converges when using it. Some pathological cases were discussed in ref. [7]. 
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4.2 Update Step for General Measures 

We now construct a general update step that converges to a local minimum of pj et (riA,Pi ) for 
a fixed set of jet constituents with momenta p^. This approach works for a wide variety of jet 
measures, including the XCone defaults. 

To motivate our general procedure, we start with the special case of the (modified) 
geometric measure, where finding a local minimum of pj e t is particularly straightforward. 
Within a given jet region A, we want to find the axis ha that minimizes 

y 'n A -Pi = n A - (y ^Pi] =n A -PA , (4.1) 

i£A i&A 

where pa = YIigaP^ is the total four-momentum of all jet constituents. Introducing a La¬ 
grange multiplier A (as in [50]), the quantity 

n A -PA + A(n^ - 1) (4.2) 

is minimized for 

n A = |^t) with p A = VVi, (4.3) 

such that the jet axis ha exactly aligns with the total three-momentum of the jet. Thus, min¬ 
imizing the modified geometric measure is equivalent to finding N mutually stable (Voronoi- 
bounded) cones. In the same way, any measure of the form 

p iet (nA,Pi) = n A - Pi fipi) (4.4) 

will be minimized by 

n A={ 1,-jfrj} with qA = ^2pif(Pi), (4.5) 

where q A is the effective total three-vector of the /-weighted jet constituents. For these cases, 
one-pass minimization will terminate in a finite number of assignment/update steps. 

The conical geometric measure does not take the form of eq. (4.4), but rather takes the 
more general form 

Pjet(nA,Pi) = n A ■ Pi g(pi, n A ), (4.6) 

where the jet measure has nonlinear dependence on n A . This means that the jet axis and the 
jet three-momentum do not in general align. For the XCone default measure in particular, 
the extra factor of cosh y A in the jet measure means that there is an offset between the axis 
and the momentum proportional to the jet mass. Thus, we cannot directly use the above 
stable-cone finding logic to minimize /0j et - Instead, as in ref. [7], we can define an update step 
based on the previous ha value: 9 

n A w = pry} with qA = ^Pi g(jH , n^ ld ). (4.7) 

9 For practical purposes, it is sometimes necessary to include an “effective mass” term by changing iia -Pi —>■ 
ua ■ Pi + e with small e to avoid potential divide-by-zero errors. 
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As long as the dependence on ua is mild enough (roughly 1 < /3 < 3 for the conical geometric 
measure), this procedure will converge within a desired accuracy in a reasonable number of 
assignment/update iterations, and we adopt this strategy for the XCone default measures. 
(In practice, due to the presence of local minima, the one-pass minimization may converge to 
a higher value of Tn than the original seed axes value. For this reason, we always return the 
smallest Tn value and associated axes seen among all update steps.) 


4.3 Seed Axes for One-Pass Minimization 


Recursive clustering algorithms are particularly effective to find seed axes for one-pass min¬ 
imization. When run in exclusive mode, a recursive clustering algorithm returns exactly X 
jets which can then be interpreted as X lightlike seed axes. In fact, the axes are often so 
good in practice that the iterative improvement step is unnecessary. One could even imagine 
a more general strategy that separates jet axis finding (here using recursive clustering) from 
jet region finding (here using X-jettiness partitions), and we plan to pursue this possibility 
in future work. Unlike generic cluster optimization, recursive clustering algorithms are com¬ 
putationally efficient, and this efficiency is inherited by our XCone implementation (at the 
expense of only guaranteeing a local Tn minimum). 

For the conical geometric measures with 7 = 1, including the XCone defaults, good seed 
axes can be found by running the generalized kx clustering algorithm with a generalized Ef 
recombination scheme. The generalized kx clustering measure [4, 63] is parametrized by an 
exponent p and a jet radius R: 

dij = min ( p%p -gf, d iB = p% (4.8) 

where p = 1 is the kx algorithm [8, 56] and p = 0 is the Cambridge/Aachen algorithm [57-59]. 
The generalized Ef recombination scheme is parametrized by an energy-weighting power 5, 
such that one obtains a massless recombined four-momentum p r given by 


PTr = VTi T PTj > 


_ Pxi&i + Pxj<t>j _ PTiVi + PxjVj 

PTi+Pxj ’ Vr P S Ti+P S Tj 


(4.9) 


where 5 = 1 is the original Ef scheme, <5 = 2 is the Ef scheme [8, 78], and <5 = oo is the 
winner-take-all scheme [27, 61, 62]. 

For finding seed axes, the recommended parameters for 0 < (3 < 2 are 




(4.10) 


with matching radius parameter R. To understand this heuristic choice, consider starting 
with a final state of N + 1 particles and running one iteration of exclusive generalized kx 
to find X axes. For this procedure to give good seed axes for Tn minimization, we want to 
choose the values of p and 5 that match the behavior of the X-jettiness metric as closely as 
possible. Essentially, we want dfB to match the beam measure Pbeam; to match the jet 
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Min. Axes for p = 2 (One pass) Min. Axes for P = 1 (One pass) 



§ 


5 


(a) 


(b) 


Figure 5: Fraction of events where all XCone jets from one-pass minimization starting from 
generalized kx jet axes as seeds align with the axes from global Tn minimization. This is 
for the BOOST 2010 top sample (Herwig 6.5, px € [500,600] GeV) [11], using the conical 
geometric measure with N = 6 and R = 0.5. (a) The XCone default (/3 = 2). (b) The recoil- 
free default (/3 = 1). Here, p and 8 parametrize the generalized kx metric and recombination 
scheme, respectively. The black boxes indicate the preferred values of p and 5 from the 
heuristic choice in eq. (4.10) (with <5 = 4 indicating 8 —» oo). 


measure pj et , and the recombination scheme to appropriately place the merged axis in the 
desired location. 

We perform this heuristic analysis for the conical measure, which is a bit easier to un¬ 
derstand than the conical geometric measure, though the same conclusions hold. To match 
the conical beam measure, generalized kx with p > 0 already gives the right behavior, since 
the softest particle farther than R from any other particle is merged with the beam. 10 

To match the conical jet measure, we want dij to depend on the combination pxiR^j, 
which is achieved for 

(4.11) 


P = 


P 


To match the conical axis behavior, we have to know which axis minimizes the Tn value for 
a jet region consisting of two particles. Labeling the two particles 1 and 2 and simplifying to 
one dimension 0 without loss of generality, we have 


Tn ~ ml0i — 4>a\ 13 + PT2 102 — 0a]' 3 , 


(4.12) 


10 In principle, it is possible to also handle the 7 ^ 1 case by further modifications of eq. (4.8) (such as those 
proposed in ref. [79]), but we have not attempted that for the present XCone implementation. 
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II 

to 

Seed axes 

One-pass min 

08 = 1) 

Seed axes 

One-pass min 

Jets 

0.95 

0.96 

Jets 

0.95 

0.97 

Events (> 4) 

0.99 

0.99 

Events (> 4) 

0.99 

0.99 

Events (> 5) 

0.92 

0.93 

Events (> 5) 

0.97 

0.98 

Events (6) 

0.78 

0.81 

Events (6) 

0.72 

0.81 


(a) 



(b) 



Table 2: Fraction of XCone jets that are aligned with the “true” minimum from global Tn 
minimization using only the seed axes from generalized kx jets and after one-pass minimiza¬ 
tion for (a) /? = 2 and (b) (3 = 1. Also shown are the fraction of all events with 4 or more, 5 
or more, and all 6 jets aligned with the global minimum. 


where 4>a is the location of the axis. Solving cYTn/^a = 0 to find the location of the 
minimum, we find 


4>a 


Pti&I +PT2<fo 

Pxi “I" PT 2 



(4.13) 


which is exactly the generalized Et recombination scheme. This is the logic behind the 
heuristic choice in eq. (4.10). 

To explicitly validate the choice in eq. (4.10), we consider a sample of boosted top quarks 
from the BOOST 2010 report [11], using N = 6 and R = 0.5. A key feature of this boosted 
top sample is the presence of initial-state radiation, which generates an additional seventh 
hard jet in the event, providing a nontrivial test scenario. We first determine by brute force 
the global Tn minimum, as best as we can, by performing one-pass minimization on a wide 
range of seed axes. Next, we perform the one-pass minimization with the generalized kx jets 
as seed axes for a range of p and 5 values. For each p and 5, we then count the fraction of 
events that have all N = 6 XCone jet axes within A R < 0.1 of the axes found from global Tn 
minimization. The results are shown in figure 5, which shows that the choice in eq. (4.10), 
shown by the black boxes, does give the best performance. We also observe that a wide range 
of 5 values give similar results, while the choice of p is more relevant, especially for (3 = 2. 

The fraction of aligned XCone jets, as well as the fraction of events where > 4, > 5, 
and all 6 XCone jets are aligned with the global minimum, both before and after one-pass 
minimization, are shown in tables 2a and 2b. Even without one-pass minimization, i.e. 
using the seed axes only, 95% of the individual jets are closely aligned with the global Tn 
minimization for both (3 = 2 and (3 = 1. This suggests that finding local Tn minima from 
generalized kx seed axes is a robust procedure that often results in a global Tn minimum. 
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The presence of additional hard jets from initial-state radiation can of course confuse 
N = 6 jet finding, leading to a roughly 70-80% success rate for correctly identifying all 6 jets 
originating from the top decays. It is also not obvious that Tn minimization will necessarily 
always yield the best boosted top reconstruction, and it might well be that “failed” Tn minima 
are still useful for physics analyses. For a detailed study of the phenomenological aspects we 
refer to ref. [1], which also explores an N = 2 x 3 strategy for this final state. 


5 IV-jettiness Factorization with Various Measures 

A key attribute that originally motivated the use of IV-jettiness is its factorization properties 
in the limit Tn —>■ 0 [5], which greatly simplifies calculations of the corresponding exclusive jet 
cross sections. The original IV-jettiness factorization theorem was derived for active-parton 
cross sections 11 using techniques from Soft-Collinear Effective Theory (SCET) [83-87], which 
we also make use of here. So far, these properties have only been fully studied for situations 
where the measure is linear in a component of the particle momenta [5, 15, 16, 21], which 
simplifies the objects appearing in the factorization theorem. The examples studied thus far 
include the geometric and geometric -R measures in table 1. 

In this section, we derive the factorization properties for more general measures. We will 
start with a generic analysis and eventually focus on /3 = 2 jet measures. We investigate the 
impact of the choice of jet axes and different beam measures. We also explain how transverse 
momentum conservation restricts the range of jet observables that can be calculated using 
the simplest version of the IV-jettiness factorization theorem. 


5.1 Separating into Jet and Beam Regions 

Due to the linear sum over particles i in eq. (2.1), IV-jettiness can be obtained by adding up 
distinct contributions from the beam and N jet regions r 

7n = = Tn + T& + Tn-\ - \~Tn ■ (5.1) 


If only a single measurement is made on the beams as in eq. (2.1), we can simply use T/y+T/y = 
T^ eam here. Thus the IV-jettiness cross section is obtained from the more fundamental cross 
section which is fully differential in the T [y for each region, 


dtr(-Wv) 

d Tn 


I n dr N s(t n ~^t 


n) 


dcr(X N ) 


d7£dT£d7^ ■ ■ • dT^ 


(5.2) 


n We only consider factorization for active-parton cross sections, initiated by incoming quarks or gluons, 
in order to avoid the complications associated with the spectator partons present for incoming hadrons, such 
as Glauber effects [80-82]. When using these active-parton factorization theorems, it is nevertheless often 
assumed that the initial-state quarks and gluons are determined by standard parton distributions. For the 
IV-jettiness observables, Glauber effects have not been fully treated in the literature. 
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where the products and sum run over r = a, b, 1,..., N. Here Xn denotes a set of measure¬ 
ments made on the N signal jets and on other final-state particles like electroweak bosons or 
nonhadronic decay products which we write as follows 


M*n) 

df^df^df^ • • • dTff 


j d<f>AT ^ S * 


da K ($N) 


dr^dr^dri 


■ ■ ■ d% N 


X N (d> 


N 


' N 


(5.3) 


Here, the sum over k runs over all relevant partonic channels k = {K a , Kb] K\, ... , kn} for the 
underlying 2 —» N process (or 2 —>• N+L where L denotes additional non-strongly-interacting 
final states). The s K is the appropriate factor to take care of symmetry factors and flavor and 
spin averaging for each partonic channel. The d<&N corresponds to the complete phase-space 
measure of the Born process with massless partons, 

fd$ N = f— — [dd> N (q a +qb]qi,-■ ■ ,qN,q)^-d$ L (q), (5.4) 

J 2 E- m J x a x b J 27r 

where dd>7v(...) on the right-hand side denotes the standard Lorentz-invariant iV-particle 
phase space, and d4>i(g) the remaining nonhadronic phase space with total momentum q. 
The variables appearing here and the restrictions we impose on the measurement function 
X N (d> N ) will be described further below. 

Now consider Tn in the exclusive iV-jet limit Tn —>• 0. Since we are interested in the 
simplest form of the factorization theorem, we assume that the jets are well separated from 
each other and from the beams, with no strong hierarchies in the jet pts. We also assume that 
if we are computing the cross section differential in 7)0, we have parametrically 7)0 ~ - 12 

For definiteness we assume that the components in the decomposition in eq. (5.5) below scale 
homogeneously, which will indeed be the case if the only iV-jettiness that we measure is the 
total Tn- In the exclusive iV-jet limit, the final state consists of only soft radiation and so- 
called n r -collinear energetic radiation which is collinear to one of the jet or beam directions 
n r . Here, the key property of iV-jettiness is the presence of the minimum in its definition, 
which leads to a linear decomposition for both Tn and Tn- Namely, they can be decomposed 
as a sum of contributions coming from each of these types of emissions, 


T n 


| '7'[ n d , q —[m] , q -fajv] T-[soft] 

' N T" ' N ' ' N “T • • • “T > N “T ' N 


.7— r _ r-[?v] I rT-r[soft] 

>N ~ >N + ' N 


(5.5) 


where the [n] superscripts refer to the contribution from emissions collinear to the n-direction, 
and [soft] to soft emissions. For definiteness, we let 


n a = (!,£), n b = (l,-z), 


(5.6) 


where z is the physical beam direction. 

Equation (5.5) encodes the fact that for all of the measures in table 1, the n r -collinear 
emissions only contribute to the measurement in the r-tli region, while the soft radiation 

12 This last assumption avoids the appearance of large nonglobal logarithms, lnfTjv/Tjv )■ These logarithms 
will not appear when considering the cross section differential only in the total Tn- 
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contributes to all regions and can itself be decomposed as in eq. (5.1). This linearity is the 
key property that allows deriving a factorization theorem which decomposes the exclusive 
IV-jet cross section into a product of functions for each type of radiation. The basic form of 
the IV-jettiness factorization theorem is [5] 

= tr 0 (g) B Kb <g> J K1 <g> ... <g> J Kn ® S% . (5.7) 

d71v 

Here, H^ is a hard function, B Ka)Kb are beam functions, J KA is a jet function for the A-tli jet 
region, and Sfo is a soft function. A description of the variables these objects depend on will 
be given below. We note immediately that Hfj depends directly on the full partonic channel 
k, as it contains the process-specific matrix elements, while S^ depends on k only via the 
color representations. The J KA depend on whether ka is a quark or gluon that initiates the 
jet, and B KatKb each depend on the flavor of the initial-state partons n a and and the type 
of initial-state hadrons. The and Sfij are both matrices in the color space of k which are 
traced over in eq. (5.7). 

The precise form of the convolutions in eq. (5.7), as well as the definitions of the beam, 
jet, and soft functions, depends on the choice of jet and beam measures used in the IV-jet tiness 
observable. On the other hand, the hard function is not affected by these choices. So far, 
we have been using the observables Tfc without specifying the method of fixing the jet axes 
n r . The form of the convolutions will generically depend on the jet axes choice. We discuss 
below the observables Tobtained after the axes minimization in eq. (2.3). The factorization 
in eq. (5.7) holds for any jet axes choice that is within 0(A) of the minimized jet axes, where 
the power counting parameter A is defined below. 

5.2 Categorizing Measures by Power Counting 

To determine the structure of the convolutions in eq. (5.7), it is first instructive to form cate¬ 
gories for the measures in table 1 that share common features in their convolution structure. 
In particular, we classify them by how they scale with the SCET power counting parameter 
A <1. Below, we use a light-cone decomposition of the momenta based on the jet axis ua 
satisfying n 2 A = 0 as well as the auxiliary vector t%a obeying n\ = 0 and ha • ua = 2. 

An n^-collinear mode within the A-th jet has momentum scaling as ( riA-Pi, nA'Pi,P^ A± ) ~ 
t%a ■ Pi (A 2 ,1, A). Here and below we use the label T to refer to components perpendicular to 
the respective jet axis Ha, while T indicates transverse momentum with respect to the beam. 
Considering all the jet measures in table 1, those with (3 = 2 have ~ A 2 (which includes 
the geometric measures), while those with (3 = 1 have ~ A. Since the components 

in the decomposition in eq. (5.5) scale homogeneously, the scaling of the corresponding soft 
momenta T)C')[ S0, ' t l must be the same as those of the corresponding collinear emissions. The 
soft momenta scale homogeneously, independent of the jet directions, so ps ~ A 2 for (3 = 2 
and Ps ^ A for [3 = 1. The (3 = 2 situation is known as an SCETj observable, while the [3 = 1 
case is referred to as an SCETn observable. 
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Since the convolutions in eq. (5.7) are always between observables with the same A- 
scaling, we can classify the jet measures by whether they are in SCETi or in SCETn. A 
similar classification can also be made for the beam measures. For collinear emissions along 
either of the two beams, 1/(2 cosh yi) ~ e — ^’1 up to power corrections. All beam measures 
having this exponential rapidity dependence are in SCETi, while those measures with just 
PTi are in SCETn. Summarizing the scaling of the measures in table 1, we have: 


SCETi jets k beams: 

SCETi jets k SCETn beams: 
SCETn jets k beams: 

SCETn jets k SCETi beams: 


Geometric (-R), Modified Geometric (-R), 
Conical Geometric (0 = 7 = 2); 

Conical (0 = 2), XCone Default; 

Conical (0 = 1), Recoil-Free Default; 

Conical Geometric (0 = 1, 7 = 2), (5.8) 


though we have not made use of the last example in this paper. 

Equation (5.5) for Tjy implies that the factorization theorem will have one convolution 
for each region it is differential in. For SCETi cases we have convolutions in (ha'P )-momenta 
between the beam/jet functions and the soft function. In contrast, for SCETn cases we 
have convolutions involving transverse or T-momenta between the beam/jet functions and 
the soft function. The homogeneous scaling for the components of A-jettiness also requires 
^ jpr'], suc h that all of the soft function convolution variables are of the same order in 
the power counting. If all jets and beams are in either SCETi or SCETn, then that theory’s 
ingredients can be used for the main components of the analysis. In the mixed case of SCETi 
jets with SCETn beams, the restriction on the radiation imposed by the measurement together 
with the power counting implies that the modes in the A-th jet can have parametrically larger 
T-momenta relative to their tia axis than the modes in the beam do relative to the beam 
axis, since p l nA± ~ A > p l nab± ~ A 2 . 

One can also derive factorization theorems for 3V-jettiness measures with generic 0. For 
any 0 such that 0 — 1 » A these measures fall in the SCETi category, and they lead to 
/3-dependent jet, beam, and soft functions. This is analogous to the factorization theorems 
derived in e + e _ —>• dijets for general angularities [88, 89] and their recoil-free variants [27]. 13 
For simplicity we will not discuss the general 0 case here, but instead focus on the represen¬ 
tative cases of 0 = 1, 2. 


5.3 Impact of Axes Minimization 

In general, the jet axes ha need not align perfectly with the jet three-momenta j/4, as long 
as the difference is 0( A). That said, the structure of the factorization theorem will simplify 
if we align the ha axes within 0( A 2 ) of the jet direction. For jets defined with XCone, 

13 In the case of recoil-free angularities, there is a smooth interpolation between SCETi and SCETn as 0 
goes from 2 to 1 [27]. 
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this alignment happens automatically for any of the /3 = 2 measures (including the XCone 
default), as explained near eq. (4.1) and discussed previously in ref. [7] (see also ref. [33, 50]). 
For this reason, we will focus the remainder of our discussion on jet measures in the SCETj 
category, including the XCone default. This minimization implies that we are now discussing 
the specific IV-jettiness observable Tn rather than the generic T/v- 

The alignment of ua with pa means that the jet momentum has 0( A 2 ) perpendicular 
momentum relative to this axis. For all the geometric jet measures the perpendicular mo¬ 
mentum is actually zero, and the component observables 7"]y then have a simple physical 
interpretation, since they measure the jet mass m 2 A for each jet region via 7(y = rri A /Q a with 
Qa = 2 pEa [19]. For our XCone default measure the perpendicular momentum is 0( A 2 ), 
however this same physical interpretation still applies, with the only difference being that 
Qa = R 2 Ea/ cosh i/A- On the other hand, for the conical geometric measure with (3 = 7 = 2 
there is not a precise relation between and m\, unless we were to adopt as an additional 
approximation yi ~ y a- 

Without aligning the jet axes and the jet three-momenta, the jet functions in the IV- 
jettiness factorization theorem would depend on both QaJ and the total p nA ±, such as 
in the jet function 

Jk a {QaT^ a] ~ Pn A ±;l- J <) ■ ( 5 - 9 ) 

Here, the two terms in J KA are both 0( A 2 ), and ka indicates a quark or gluon. With the 
axes minimization, the dependence on the transverse momentum drops out, and this becomes 
simply 

J KA (Q A Tp A] ,p). (5.10) 

These jet functions, which appear in the IV-jettiness factorization theorem, are inclusive 
because the collinear radiation is always completely contained in the corresponding jet region. 
This means that they are a function of a single variable and do not depend on the jet boundary. 
However, the type of inclusive jet function we have does still depend on the jet measure. For 
instance, the geometric measures yield the standard inclusive hemisphere jet function, but 
we obtain a different inclusive jet function for the /3 = 7 = 2 conical geometric measure. 

5.4 Hadronic and Partonic Momentum Conservation 

The remaining ingredients that influence the form of the factorization theorem are momentum 
conservation and the choice of measurements Xjy made on the jets and the nonhadronic 
particles. We will discuss the first issue here, before explaining why they impact the structure 
of the factorization theorem in the next subsection. 

Momentum conservation says that 

Pbeam = Pa + Pb = ^ ’ t 5 ' 11 ) 

A 

where p A is the sum of all four-momenta for particles in region A. the p£ b (bc am ) include the 
incoming proton momentum (momenta) minus the sum of the outgoing momentum of particles 
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in the associated beam region, and q^ is the total outgoing momentum of any nonhadronic 
particles. Even if the N -jettiness measurement specihes only a single beam region, we can 
divide the beam region in two by making an artificial split at zero rapidity into regions a and 
b. This split is useful for the discussion below, since it makes it simpler to talk about the two 
beam functions that are important for the dynamics of the beam region. We set q^ 1 = 0 for 
cases where the final state does not involve nonhadronic particles. 

The largest 0(A°) momentum component from each jet and beam region in eq. (5.11) 
can be extracted by projecting along the associated IV-jettiness axis, 

11 
71 ' 

# = Wr ^ + 0( A). (5.12) 

This determines the variables appearing in the hard function 

H N = H N ({uj r n r },q,ij), (5.13) 

where r runs over a, b, 1, ..., N in the set of variables in {•••}. 14 These phase-space 
variables include things like the transverse momentum and rapidity t]a of each jet, as 
well as the overall rapidity of all non-forward radiation Y which determines the boost of the 
partonic hard collision relative to the center-of-mass frame. These hard-function variables 
form the basis for the measurements we make on the jets as specified by Wv(d > iv) in eq. (5.3) 
where q r = u r n r /2. The variables are not all independent, since momentum conservation 
correlates the large 0(A°) components of eq. (5.11). This is the same as imposing momentum 
conservation for the underlying hard partonic process with incoming and outgoing massless 
partons, 

n LL u 

Tin 'Flu n X —^1 A 

Wa Y +Ub it =q ‘ ( 5 ' 14 ) 

A 

In particular, this formula is used to compute H^r when integrating out hard modes by 
matching QCD to SCET using calculations of S-matrix elements in the two theories. And 
this momentum conservation appears above in d<hjv hi eq. (5.4). The same hard function in 
eq. (5.13) appears in the factorization theorem for exclusive jet cross sections for all choices 
of the JV-jettiness jet and beam measures. 

In eq. (5.14), the beam variables can be rewritten in terms of the total center-of-mass 
energy E cm and momentum fractions x a b for the colliding partons in the hard collision via 
uj a = x a E cm and ojb = XbE cm . The jet variables oja are chosen so that oja = 2 Ea + 0(A), 
where Ea is the true jet energy, and the presence of 0(A) contributions in this relation 
ensure that eq. (5.14) is exactly satisfied. The presence of these 0(A) terms does not affect 
the evaluation of the hard function in eq. (5.13), where we may simply replace oja —> 2 Ea- 
This same replacement should be made in the formulas for the Qa factors appearing in the 

14 To emphasize that Hn can always be written in terms of Lorentz-invariant phase-space variables, one can 
rewrite this as -ffjv({oi r uv n r ■ rv}, {nv n r ■ q}, q 2 ) with r and r' running over a, b, 1, ..., N. 
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(Modified) Geometric 

Geometric -R 

Modified Geometric -R 

XCone Default 

Qa = Po ua 

p{R, va) va 

Pc{R,Va)u a 

R 2 

- <* h 4 

2 cosh yA 


Table 3: Values of Qa for various measures. The approximation oja = 2Ea is valid as long 
as the same replacement is made in the hard function. 


jet functions, which are otherwise given by the results in table 3. However, the 0(A) terms 
can have implications for the convolutions between the jet, beam, and soft functions. To see 
explicitly how these 0(A) terms arise, it is convenient to project eq. (5.14) both along and 
transverse to the beam axis, giving 

u a = n b - q + 2 _^uja —-— 

A 

0 = 2 q% + ^ oja n AT . 

A 

The two equalities in eq. (5.15) simply fix oj a , b regardless of how precisely we specify the jet 
axes ha, the jet variables oja, or q^. This leaves the two constraints from eq. (5.16), which will 
be very important in the next subsection. These constraints involve n AT , which is determined 
by the azimuthal angle cj) nA for the axis of each jet region, but they do not depend on the 
longitudinal (rapidity) component of ha- 

5.5 Convolutions from Transverse Momentum Recoil 

We now show how the two constraints in eq. (5.16) can influence the form of the convolutions 
appearing in the factorization theorem. Throughout this discussion, we assume that the jet 
axes riA and jet three-momenta pa are perfectly aligned, as is the case for the /3 = 2 measures 
with the minimized T/v- We start with pure SCETj observables before mentioning what 
happens with SCETn beam measures. 

To begin, imagine making highly granular measurements of the jet energies and directions 
with very fine p^, Pa, and (j)A bins, as well as fully measuring the nonhadronic q ft . In this 
situation, we have effectively completely measured the transverse vector q£, the jet energies 
Ea, and the vectors tia, so we actually have a measurement that is sensitive to the 0(A) 
amount by which the oja variables differ from 2 Ea- Here, the H-th jet’s momentum can be 
written as 

Pa = {2E A -n-p A )^Y +n-p A ^Y + P A j_, (5.17) 

where n ■ pa ~ Tjy ~ A 2 . We can therefore see that the components beyond Eati a are 
0(A 2 ) and do not have 0(A) projections on the axis transverse to the beam. If we consider 


v-^ n a • riA 

oj b = n a - q + 2 _^oj a --—, (5.15) 

A 

(5.16) 
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transverse momentum conservation using the original momentum conservation in eq. (5.11), 
and insert eq. (5.17), then we find that the balance of transverse momenta at 0(A) is given 
by 

k T = PaT +PbT = + Ea U AT • ( 5 - 18 ) 

A 

Using eq. (5.16) we can see that this is a small momentum kt^ ~ A. For the beam variables p^ b , 
these 0(A) transverse components come from the transverse momenta of radiation emitted 
in the beam regions (since the transverse momenta in the proton are ~ Aqcd which is much 
smaller). For the jet components, this 0(A) momentum comes from the mismatch between 
loa and 2Ea, which we can see explicitly by using eq. (5.17) in eq. (5.16) to give 

k r = Yl { Ea - ir) n ^ ■ ( 5 - 19 ) 

A 

With the assumptions above, the constraint in eq. (5.18) is present because by making 
such a granular measurement, we have indirectly measured kj ,, and hence the total transverse 
momentum recoil of the beam radiation. This measurement therefore leads to p^-dependent 
beam functions in the factorization theorem, which appear as 

J d 2 p T B Ka (t a ,x a ,PT, B Kh (t b ,x b , k T -p T ,p) ■ (5.20) 

Here t a = oj a T^ la ^ and t b = uj b T^ involve the variables that are convolved with the soft func¬ 
tion. The double differential beam functions B Ka (t a ,x a ,px, p) were discussed in refs. [90, 91]. 
In ref. [21], examples where transverse momentum convolutions connect a jet and beam func¬ 
tion were discussed for an SCETi type 1-jettiness in deep inelastic scattering, and eq. (5.20) 
is the analog of the center-of-mass 1-jettiness variable considered there, except with the jet 
function replaced by a second beam function. The double differential factorization theorem 
with an explicit measurement of 0-jettiness and kx in SCETj was derived in ref. [92], and 
involves precisely the combination in eq. (5.20). 

To obtain a simpler factorization theorem that does not involve p-r-dependent beam 
functions, we just have to perform a less granular measurement that does not constrain every 
aspect of the final state. For cases with external nonhadronic particles, the simplest approach 
is to not fully constrain all components of q£, for example by specifying qx only within a 
bin centered on g|? ntral with width > Ag“ ntral . Since A ~ tua/Ea — 0.1, this corresponds to 
the typical size of bins that are already used in experimental analyses (unless they are only 
interested in measuring qx). This method was used in ref. [15] when deriving the active-parton 
factorization theorem for beam thrust or 0-jettiness, where qx was simply not measured. For 
beam thrust there are no jets, so g^ ntral = 0, but this approach works equally well for (N > 1)- 
jettiness where g“ ntral ~ A 0 is large. Once one uses this coarser qx binning, there are no other 
0(A) constraints on the transverse momenta. In particular, specifying the bin for qx yields 
an additional unrestricted integration over kt^ which appears in eq. (5.20) when deriving 
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the factorization theorem. Therefore, we obtain independent transverse integrals over the 
two beam functions, f d 2 pxB K (t, x,px, p) = B K (t,x, p), and only these /^-independent beam 
functions appear in the A r -jettiness factorization theorem, as in 

Bk b {tai *E<xj M) B Kb (tbi Xfj, p) • (5.21) 

Alternatively, for cases where N > 2, we can exploit the fact that we do not need to 
make finely-binned measurements of the jet energies or jet pts. We can instead be satisfied 
with a measurement with center p“ ntral in a bin of width > A/>S]? ntral , which could be for 
example using a bin centered at 500 GeV with width 50 GeV. This can be applied to both 
cases with (q ^ 0) or without (q = 0) additional nonhadronic particles. Since we can now 
vary by 0(A) at least two of the oja variables, we again loosen the constraint fixing kj- and we 
can again freely integrate over this variable, and hence also obtain eq. (5.21). Both of these 
approaches to obtaining the simpler form of beam functions in eq. (5.21) require making less 
granular measurements when specifying Xjy, but still remain fully sufficient for all standard 
LHC jet-style measurements. The only cases where eq. (5.20) become relevant is if we are 
actually interested in making a jet measurement so finely-binned that we can infer the small 
PT spectrum of the beam radiation. 

Just like for the jet function in eq. (5.10), the beam functions in eq. (5.21) are inclusive 
because collinear radiation along the beam directions is completely contained in the beam 
regions. Thus, they do not depend on the boundaries between the beam and jet regions. 
In principle, they could still depend on the beam measure, but because of eq. (3.11), for all 
the SCETi beam measures we consider here, they are always given by the standard inclusive 
hemisphere beam functions [15, 93]. 

It is interesting to consider how the above arguments change if we maintain SCETj 
measures for the jets (and aligned jet axes obtained from minimization) but now consider a 
SCETn measure for the beam; this is the case encountered in the XCone default measure. In 
this situation, we still have inclusive jet functions that do not depend on p nA ± as in eq. (5.10). 
The key change is that now the IV-jettiness measurement forces the beam transverse momenta 
to be smaller, p^ T ~ pk T ~ A 2 , and the resulting SCETn beam functions are of the broadening 
variety with t a = and tb = 7variables that are themselves 0{ A 2 ). In addition to the 
renormalization scale p, the beam functions depend on a rapidity renormalization scale v, in 
the combination v/oj a b- The v scale is needed to sum logarithms associated with rapidity 
divergences that appear from the separation of modes in the beam and soft functions [94, 95]. 
From eq. (5.18), we must also have kj, = p^ T + p^ T = q£ + YIa Ba^at ~ ^ 2 - O nce again we 
can integrate over ki £ either by considering a bin for or a bin for two of the jet energies 
Ea ■ In this case, the bins need only have a size of > A 2 g“ ntral or > A 2 £A) entral in order to 
sufficiently integrate over fcj. such that we get /^-independent beam functions, as in eq. (5.21). 

5.6 Factorization Theorems for A-jettiness 

We now have all the ingredients needed to assemble the factorization theorem for IV-jettiness 
for various jet and beam measures. For jet and beam measures in SCETi, the mathematical 


- 29 


derivation of this factorization theorem follows closely the detailed derivation given for beam 
thrust in ref. [15], or for DIS 1-jettiness in ref. [21], which we therefore will not bother to 
repeat here. The IV-jet case has also been discussed in some detail in refs. [16, 30]. The 
required ingredients in the derivation have all been discussed in the previous subsections. 

With jet and beam measures in the SCETi category, axes determined by minimization, 
and the choice for Xjy that does not directly or indirectly measure the transverse momentum 
of the beam radiation, the factorization theorem in eq. (5.7) becomes 


da K («hjv) 


dTfi&T'&T}, ■ ■ ■ d7f 


(5.22) 


= tr H^({u r n r }, q, /j,) f [jj dT^ r] u a B Ka (u a T^ a] ,x a , fi) uj b B Kb (u b T^ b] , x b , n) 

** r 

X Qi J K1 {QxT [ ; il] ,»)■■■ Q N J KN ( Qn T$ n] , /i) S% ( {U - rf? r] } , {^ }, fi) , 


where r and r' = a, 6,1,..., A" and all of the convolutions are now made explicit. Here, 
the soft function S 1 ^ depends on the IV + 2 observables 7)0. It is a scalar function of the 
variables {uj r rir /Q r }, which encode the dependence on the angles between various beam and 
jet directions through their dot products. Although not indicated by our notation, the soft 
function also depends on the size and shape of the jet regions through the precise definition of 
the jet and beam measures used to define these observables. Both the jet functions and beam 
functions in eq. (5.22) are of the inclusive variety, and hence do not depend on the boundaries 
between the jet or beam regions. The beam functions also contain the nonperturbative parton 
distributions _/)•(£, fj,) through a factorization from the perturbative radiation into calculable 
coefficients X l3 [15, 93, 96], 


_^ r 1 \ 

Bi{uk,x,n) = 

J X Z \ Z / 


(5.23) 


With geometric (and related) measures, eq. (5.22) was the version of the 1-jettiness factor¬ 
ization theorem used for the analysis in Ref. [19]. 

For the various geometric measures, the Qa factors needed for eq. (5.22) are given above 
in table 3. For the j3 = 7 = 2 conical geometric measure we let Qa = R 2 ^a- For this 
measure, the inclusive jet functions become Jk, a {QaTn 1a \ ua, h) e T (5-25), due to the 
cosh da/ coshyj weighting factor in the jet measure. Thus, they are not just the standard 
hemisphere jet functions. Similarly, for this case we also will have a soft function that can 
depend on the ua variables. 

In eq. (5.22) we are differential in two beam regions, 7jy and 7jy. If we only want to 
consider a single beam region and measurement observable 7)y eam = + 7)^, then it is 

possible to simplify the form of the factorization theorem. Using the corresponding collinear 
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projection, = 7/ beam ^, yields a “double-beam function” for SCETi measures 

BBij ( u a T ^ beaml , Wfo T^ nbeam] , x a , x b , p) = uj a uj b j d k B { (u a k, x a , p) Bj (cu fe (7/ beaml - k), x b , p) . 

(5.24) 


• , • ,i r, r ,• • ,i • /-i-afsoftl . /-^[softl ,-,-beam [soft 1 ,, • i 

Projecting the sort function m the same way, using / N L J + / N L 1 = I N , this reduces 

eq. (5.22) to 


dcr K (4>7v) 


dT/ am dT/--dT/ 


(5.25) 


= tr Hfr({uj r n r },q,n) / [[pT/ 


BB, 


K'dK’b 


(o; a ry beaml , OJ h T [ ^ m] ,X a ,X b ,p) 


X QiJ Kl (Qi T/ l] ,//)••• Qat J KJV (Q;v rjr 1 , /r) S# I} [{U - T [ ^ ] }, {^}, m) 


where now r = beam, 1,... ,1V. For the modified geometric(-R) measure, the soft function 
Sp"* in eq. (5.25) has a C-parameter-type measurement for its 7j\) eam observable and thrust- 
type measurements for the jet observables 7/, and eq. (5.25) involves the standard inclusive 
hemisphere jet functions. 

Next, we consider the mixed measure case, with SCETi jet measures and SCETn beam 
measures, still with jet axes determined by minimization and a choice of Xn that is insen¬ 
sitive to transverse momentum of the beam radiation. For this case, there has not yet been 
any literature providing a detailed mathematical derivation of a factorization theorem. Fac¬ 
torization theorems have been worked out for pure SCETn measurements of event shapes in 
e + e - —>• dijets [27, 94, 95, 97, 98], and active-parton factorization theorems have also been 
derived for pp —)• H with an Et jet veto [99] or p!^ veto [99-103]; see also [104] for transverse 
thrust. Experience from these results enables us to anticipate the form of the convolutions 
that will appear between the beam and soft functions in the mixed measure IV-jettiness case. 
So even though the complete derivation of the factorization theorem for this case is beyond 
the scope of this work, we can still put the information collected above together to anticipate 
its structure. 

For SCETn beam measures we expect the double-beam function to be given by 


BBij( k, x a ,x b , p,—, — ) = /d k' Bi(k',x a ,fJ,,—) Bjlk- k',Xb,n, — ) ■ (5.26) 

UJ a U} b J J V u a J J \ UJ b J 

The individual beam functions here are of the broadening type and involve the rapidity scale 
parameter v [95]. For SCETi jet measures and a single SCETn beam measure we then expect 
a factorization theorem of the form 


du K (<I> 


N) 


cir^dri • • • dr/ 


= tr H^({u r n r },q,p) j\j^ d7j! 


n r \ 

N 


bb Kc 


Kb 


/ T~L^'beamJ 

7 N 


V V 


%a5 M? 


UJ a UJ b 


(5.27) 
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X Qi j k1 (Qi rir i] , /i) • • • QtvJ KJV (Qtv r^ 1 , /i) s# I/n) ({7^ - rjT r] }, {^^. 

This is the factorization formula that is relevant for the XCone default measure, with the 
Qa factors given above in table 3. Note that here the soft function 5^ I//n ^ has broadening- 
type variables convolved with the beam functions, and has dependence on the scale u which 
compensates the v dependence in the double-beam function. The conical measure with /3 = 2 
will have an analogous factorization theorem but requires different jet and soft functions that 
take into account that the jet measure cannot be written as n- Piuja/Qa with some Qa- We 
leave a detailed mathematical analysis and proof of the active-parton factorization theorem 
in eq. (5.27) to future work. It will also be interesting to test it against fixed-order predictions 
for these IV-jettiness distributions. 

The other main class of measures in table 1 are those that have both jet and beam 
measures in the SCETn category. This includes the recoil-free default XCone measure, as well 
as the conical measure with (3 = 1. Once again there has not yet been a detailed mathematical 
analysis of this case in the literature, but from our previous analysis and from experience with 
simpler cases, we can anticipate the form of the associated factorization theorem. With axes 
determined by minimization, and with a choice of Xn that is again insensitive to the total 
transverse momentum of the beam radiation, we expect the appropriate factorization theorem 
to contain the same SCETn double-beam function in eq. (5.26) with no additional recoil 
convolutions. This would be analogous to the factorization theorem for recoil-free broadening 
in e + e _ collisions in ref. [27]. In contrast to eq. (5.27), the jet functions must now be of the 
broadening type and likely also depend on a rapidity scale u. The corresponding soft function 
Sf l) now only depends on convolution variables of the broadening type and has to cancel 
the v dependence of both beam and jet functions. We again leave a detailed mathematical 
analysis and proof of the active-parton factorization theorem for this case to future work. 

6 Conclusions 

In this paper, we introduced the new XCone jet algorithm, which is based on the IV-jettiness 
event shape. XCone is an exclusive cone algorithm that finds a fixed predefined number 
of jets. Exploiting the measure flexibility inherent to IV-jettiness, we defined a new conical 
geometric measure that combines the geometric measure, which is theoretically motivated 
and preferred, with the conical measure, which has already been proven to be experimentally 
robust in the context of jet substructure techniques using IV-subjettiness. In a companion 
paper [1], we present three physics case studies to highlight how XCone can be beneficial 
to a variety of LHC analyses. In particular, XCone is capable of resolving overlapping jets 
without requiring a separate split/merge step, and allows for a continuous transition from 
the resolved regime of well separated jets to the boosted regime of overlapping jets. 

Our focus in this paper was on the case 7 = 1, for which the beam measure scales as 
pr, such that T/v minimization is roughly the same as minimizing the total unclustered px- 
By changing 7 , one changes whether jets are found preferentially in the central or forward 
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parts of the detector. In the future, it would be interesting to study the impact and utility 
of different 7 values, especially 7 = 2 which is the natural value from the original geometric 
measure. At present, the XCone code is limited to 7 = 1, primarily because our method to 
find seed axes employs the existing longitudinally-invariant generalized kx algorithm. It is 
possible to build recursive clustering algorithms optimized to find seed axes for any given Tn 
measure, which is planned for future work. 

In constructing the XCone algorithm, we have chosen a specific measure for both the N- 
jettiness partitioning into jet and beam regions as well as the jet axis finding via the overall 
IV-jettiness minimization. This has lead to an interesting compromise, where in order for the 
XCone default measure to use dot-product distances in the jet partitioning, the jet regions 
could not be perfectly stable cones (meaning the jet axis is not exactly aligned with the total 
jet momentum). One could imagine loosening the requirement of Tn minimization, though, 
to define an array of exclusive jet algorithms. Following the idea that jet axis finding and jet 
region finding can be regarded as two distinct steps, one could use any exclusive clustering 
algorithm to find jet axes and only use Tn for defining the jet partitions. Alternatively, if 
one wants the jet axis to be perfectly aligned with the jet momentum, one could build an 
exclusive cone jet algorithm that directly searches for N mutually stable perfect cones. More 
generally, it is worth reexamining the potential of exclusive jet algorithms at hadron colliders, 
and XCone provides a clear proof of concept with interesting physics applications [1]. 

Beyond just being an exclusive jet algorithm that finds a fixed number of jets, XCone 
can be adapted to become an inclusive jet algorithm that finds a variable number of jets by 
analyzing the distribution of Tn for different N. For an event with M jets, Tn should be 
large when N < M and small when N > M, producing a sharp downward transition in the 
value of Tn when N = M. Therefore, one could iteratively increase the value of N until Tn 
undergoes this transition, either by measuring the “slope” dT/v/dlV or by imposing a fixed 
Tut- 15 Using XCone as an inclusive jet algorithm could potentially be useful for jet counting 
in event samples with a variable number of jets, for accurate event reconstruction in the face 
of hard initial state radiation, or for improving background discrimination by dividing an 
event sample into exclusive IV-jet bins. 

Finally, we anticipate that the XCone default measure will be used in future IV-jettiness 
theoretical calculations. Since XCone is IRC safe, there are no obstacles for performing 
fixed-order or resummed calculations for any of the measures studied here. While jet and 
beam measures that are linear in the particle momenta (like the XCone default measure) are 
simplest when using factorization to carry out calculations, the discussion in section 5 implies 
that the same SCET-based methods can also be applied for other measures. Ultimately, 
we look forward to comparing precision XCone-based calculations to precision XCone-based 
measurements at the LHC. 

15 Because XCone only finds a local minimum by default, there is no guarantee that Tn is a strictly decreasing 
function of N, though in practice this is a small effect when using the heuristic in eq. (4.10). 
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